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Abstract 

A high statistics sample of photoproduced charm particles from the FOCUS (E831) 
experiment at Fermilab has been used to measure the D° and D + lifetimes. Using 
about 210 000 D° and 110 000 D + events we obtained the following values: 409.6 ± 
1.1 (statistical) ± 1.5 (systematic) fs for D° and 1039.4 ±4.3 (statistical) ± 7.0 (sys- 
tematic) fs for D + . 



The study of the charm hadron lifetimes has been fundamental for our under- 
standing of the heavy quark decays. The most important contribution is the 
spectator diagram which contributes equally to the widths of all hadrons of a 
given flavour [1]. In the early days of charm physics it was quite a surprise when 
the experiments measured a large value for / td° ■ It is generally believed 
that this large ratio (~ 2.5) is mainly due to the destructive interferences 
between different quark diagrams that contribute only to D + decays. The in- 
creasingly precise measurements of the heavy quark lifetimes have stimulated 
the development of theoretical models, like the Heavy Quark Theory [2] , which 
are able to predict successfully the rich pattern of charm hadron lifetimes, that 
span one order of magnitude from the longest lived (D + ) to the shortest lived 
(fi°). 

In this letter we present the most accurate measurement to date of the lifetimes 
of the D + and D°. Although the accuracy reached by the previous experiments 
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is remarkable, for example the D° lifetime is known with an uncertainty of 
~ 1% [3], we think that a more precise determination of the D° lifetime 
would be needed. For example it would allow a more accurate check for the 
determination of the lifetime difference in the neutral D-meson system (to 
evaluate the parameter y = Ar/2r of the D° — D° mixing [4]). 

Charmed particles were produced by the interaction of high energy photons, 
obtained by means of bremsstrahlung of electron and positron beams (with 
typically 300 GeV endpoint energy), with a beryllium oxide target. The mean 
energy of the photon beam was approximately 180 GeV. The data were col- 
lected at Fermilab during the 1996-97 fixed-target run. More than 6.3 x 10 9 
triggers were collected from which more than 1 million charmed particles have 
been reconstructed. 

The particles from the interaction are detected in a large-aperture magnetic 
spectrometer with excellent vertex measurement, particle identification and 
calorimetric capabilities. The vertex detector consists of two systems of silicon 
microvertex detectors. The upstream system consists of 4 planes interleaved 
with the experimental target [5] (2 target slab upstream then 2 silicon planes 
and the replica of this setting), while the downstream system consists of 12 
planes of microstrips arranged in three views. These detectors provide high 
resolution separation of primary (production) and secondary (decay) vertices 
with an average proper time resolution of ~ 35 fs. The momentum of the 
charged particles is determined by measuring their deflections in two analy- 
sis magnets of opposite polarity with five stations of multiwire proportional 
chambers. Kaons and pions in the D-meson final states are well separated up 
to 60 GeV/c of momentum using three multicell threshold Cerenkov counters. 

The final states are selected using a candidate driven vertex algorithm [6]. A 
secondary vertex is formed from the reconstructed tracks and the momentum 
vector of the D candidate is used as a seed to intersect the other tracks in 
the event to find the primary vertex. Once the production and decay vertices 
are determined, the distance i between them and the relative error ae are 
computed. Cuts on the E/at ratio are applied to extract the D signals from 
the prompt background. The primary and secondary vertex are required to 
have a confidence level greater than 1%. 

The vertices (primary and secondary) have to lie inside a fiducial volume^ 
and the primary vertex must be formed with at least two reconstructed tracks 
in addition to the seed track. The Cerenkov particle identification cuts used 
in FOCUS are based on likelihood ratios between the various stable particle 



The reason for this cut is the presence of a trigger counter just upstream of 
the second microstrip device, therefore we define the fiducial volume as the region 
between the first slab of the experimental target and this trigger counter. 
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identification hypotheses. These likelihoods are computed for a given track 
from the observed firing response (on or off) of all cells within the track's 
(f3 = 1) Cerenkov cone for each of our three Cerenkov counters. The product 
of all firing probabilities for all cells within the three Cerenkov cones produces 
a x 2 -hke variable Wi = — 2 ln(Likelihood) where % ranges over the electron, 
pion, kaon and proton hypotheses (see Ref. [7] for more details). We require 
Ak = W n — Wk > 1, called kaonicity, for the tracks reconstructed as a 
kaon. Analogously the tracks reconstructed as pions have a pionicity, A w = 
W K — W w , exceeding I. 

In Fig. la and lb we show the invariant mass plots obtained with this set 
of cuts and with tjot > 9 for the decay modes D° — > K~n + and D° — > 
K~ty + 'k + 'k~ respectively (throughout this paper the charge conjugate state is 
implied). Fig. lc is the invariant mass plot for the decay mode D + — > K~-n + -n + 
with tjoi > 14. This set of cuts is chosen to optimize the yield and the 
background underneath the signal (S/N ratio). 

From a binned maximum likelihood fit we find 139 433 ± 520 D° — > K~tt + , 
68 274±360 D° -> ^-tt+tt+tt- and 109 877±385 D + -> K-ti + tv + candidates. 
The plots are fit with two Gaussians with the same mean but different widths 
to take into account the different resolution in momentum of the tracks passing 
through one or two magnets (see Ref. [6] for more details) of our spectrometer 

plus a 2 order polynomial. The low mass region is excluded in the fit to 
avoid possible contamination due to other hadronic charm decays involving 
an additional n°. 

The lifetime is measured using a binned maximum likelihood fitting tech- 
nique [8]. A fit is made to the reduced proper time distribution in the signal 
region. The reduced proper time is defined by t' — (£ — Nag)/ (fijc) where £ 
is the distance between the primary and the secondary vertex, ag is the reso- 
lution on £ and N is the minimum "detachment" cut required to extract the 
signal. Our vertexing algorithm provides very uniform reduced proper time 
acceptance even at very low reduced proper times. If absorption and accep- 
tance corrections are small enough that they can be neglected, and if ag is 
independent of £, one can show that the t' distribution for decaying charmed 
particles will follow an exponential distribution. These assumptions are very 
nearly true in FOCUS [6]. 

The signal region reduced proper time distributions (indicated with dashed 
lines in Fig. 1) are formed from events with invariant mass within ±2cf of the 
mean D mass; a stands for weighted sigma because the invariant mass plots 
were fitted with two Gaussians. The dependence of the lifetime measurement 
on the choice of the signal and background region is discussed later in the text. 
The binned maximum likelihood method allows direct use of the proper time 
distribution of the data above and below the D mass peak to represent the 



4 



background underneath the signal instead of using a background parametriza- 
tion. We have chosen two sidebands starting 4cf above and below the mean 
D mass, each half as wide as the signal region (indicated with dotted lines 
in Fig. 1). The signal and background reduced proper time distributions are 
binned in proper time wide bins (200 fs) spanning about 10 nominal lifetimes. 

The observed numbers of events in a reduced proper time bin % (centered at 
£•) in the signal and sideband histograms are labeled and bi respectively. 
The predicted number of events rii in a reduced proper time bin is given by: 

n>i - (iV s - B) — —— T - + B—— (L) 

where N s is the total number of events in the signal region, B is the total num- 
ber of background events in the signal region and f{t'j) is a correction function. 
The fit parameters are B and r. The /(£•) correction function, derived from 
a Monte Carlo simulation, corrects the reduced proper time evolution of the 
signal for the effects of geometric acceptance, reconstruction efficiency, anal- 
ysis cuts, hadronic absorption and decay of charm secondaries. The use of a 
multiplicative /(^) correction, rather than an integral over a resolution factor, 
is justified since our reduced proper time resolution (~ 35 fs) is much less than 
the D° or D + lifetime. 

A separate /(£') correction function is used for each of the three decay modes. 
Our Monte Carlo simulation includes the Pythia [9] model for photon-gluon 
fusion and incorporates a complete simulation at the digitization level of all 
detector and trigger systems and includes all known multiple scattering and 
particle absorption effects. The Monte Carlo was run with ~ 15 times the 
statistics of the experiment. 

The plots a) of Fig. 2, 3 and 4 show the correction function f(t') for the three 
decay modes in bins of reduced proper time. The f(t') function is obtained 
by dividing the simulated reconstructed charm yield in each bin by the input 
decay exponential. The fall off in f(t') for the D + case is due to the exclusion 
of long lived events with vertices downstream of the fiducial volume 1 . 

A factor is included in the likelihood function in order to relate B to 
the number of background events expected from the side band population. 
The background level is thereby jointly determined from the invariant mass 
distribution and from the reduced proper time evolution in the side bands. 
The likelihood function is then given by: 

^ = ^signal x ^bg ( 2 ) 



5 



where 



_ b jy ^exp(-n t ) 
-signal — 11 | W 



and 



4>g = 



5 b g 



bg ! 



exp(— £>) 



(4) 



with -A^bg = Si &i ( we assume a linear background because the 2 order term 
of the polynomial is negligible). 

The plots b) of Fig. 2, 3 and 4 show the predicted events (histogram) superim- 
posed on the observed events, the background events hi are also superimposed. 
In plots c) of Fig. 2, 3 and 4 a pure exponential function with the fitted life- 
time is superimposed on the background subtracted and f(t') corrected t' 
distribution. 

The measured lifetimes are: 408.75 ± 1.42 fs for D° -> K~ir + , 411.25 ± 1.95 fs 
for D° -> K-iT+iT+iT- and 1039.42 ± 4.28 fs for D + -> K-it + it + . 

Our lifetime measurements have been tested by modifying each of the vertex 
and Cerenkov cuts individually. For example in Fig. 5 one can see the measured 
lifetimes versus the £/ag detachment cut. The measured lifetimes of the three 
decay modes are stable with respect to £/&£. Our measured lifetimes show no 
significant variation with the cuts employed to extract the signal. 

We check our D° lifetime evaluation partitioning the total sample into D*— tag 
and no-tag according to their origin. The obtained lifetimes are in very good 
agreement with the reported values. 

To further check our lifetime measurements we have used tight cuts in order to 
extract a signal with virtually no background. Fig. 6 shows the invariant mass 
plots of the three decay modes for this set of cuts. The lifetime measurements 
from these samples, 411.26 ±3.11 fs for D° -> K~n + 413.10 ± 4.80 fs for 
D° -> R-n+Ti+Ti- and 1036.66 ±8.00 fs for D+ -> are in very good 

agreement with our previous determination. 

Systematic uncertainties for these lifetime measurements can arise from several 
sources. We performed a detailed study to analyze these sources. 

There is an uncertainty due to the absolute time scale which was determined 
by studying the absolute length and momentum scale in the experiment (see 
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Ref. [10] for more details). We estimate an uncertainty of ±0.11% for this 
source. 

Another source of systematic uncertainty is linked to the detector and recon- 
struction efficiency The f(t') corrects the reduced proper time distribution 
for these effects, but an uncertainty could originate if there is a mismatch 
between the Monte Carlo simulation and the data. We have verified that our 
Monte Carlo accurately reproduces the distributions of several relevant vari- 
ables, such as the longitudinal and transverse momenta, the multiplicity of 
the production vertex, the measured decay length and the estimated error 
on the reconstructed proper time. In order to estimate this uncertainty we 
split our total sample into independent subsamples depending on D momen- 
tum, particle versus anti-particle and the different periods in which the data 
were collected. The splits into D momentum and charge conjugation are the 
natural tests to reveal a possible mismatch between data and Monte Carlo be- 
cause they probe the response of the detector. The main reason for the period 
dependence is the insertion of the upstream silicon system (which improved 
the resolution) in the target region during the 1997 fixed-target run period. 
A technique, employed in FOCUS and in the predecessor experiment E687, 
modeled after the S-f actor method from the Particle Data Group [3], was used 
to try to separate true systematic variations from statistical fluctuations. The 
lifetime is evaluated for each of the 8 (= 2 3 ) statistically independent subsam- 
ples and a scaled variance is calculated; the split sample variance is defined as 
the difference between the reported statistical variance and the scaled variance 
if the scaled variance exceeds the statistical variance. This contribution to the 
systematic error is reported as split sample in Table 1. 

The reported lifetimes are obtained with a particular set of fitting conditions. 
For example the width of the bins or the range of the if distribution. This is a 
particular choice and the lifetime should be independent of it. We investigated 
if this could be a possible source of uncertainty by varying the width of the 
bins, the upper limit of the t' distribution, the location and the width of the 
sidebands, and the width of the signal region. In addition we studied the effect 
of using only the low or only the high mass sideband as well as the effect of 
eliminating the background term in the likelihood (the second term in equation 
2). For all these fit variants the sample variance is used as an estimate of this 
uncertainty because the various measurements are all taken as a priori likely. 

A further source of systematic error can be due to uncertainties in the target 
absorption corrections. Two effects are present: hadronic absorption of decay 
daughters which would increase the fitted lifetime if neglected and absorption 
of the D in the target which would tend to decrease the fitted lifetime if 
not taken into account. In our Monte Carlo all known particle absorption 
effects (values from Particle Data Group [3]) have been simulated for the decay 
daughters. For the charm hadron absorption our simulation assumes 1/2 of the 
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Table 1 



Contributions in percent to the systematic uncertainty. 





D° 


D° 


D° 


D+ 


Source 


K-TT+ 




combined 


K-TT+TT+ 


Absolute time scale 


0.11% 


0.11% 


0.11% 


0.11% 


Split sample 


0. 


0.37% 


0.13% 


0. 


Fit variant 


0.19% 


0.14% 


0.17% 


0.15% 


Absorption 


0.11% 


0.20% 


0.20% 


0.38% 


Acceptance 


0.21% 


0.21% 


0.21% 


0.52% 


Total systematic error 


0.32% 


0.50% 


0.38% 


0.67% 



cross section for neutrons. We estimate the uncertainty of this contribution by 
varying the charm cross-section by 50% and the daughter particle interaction 
cross-sections by 25% in the Monte Carlo simulation. We verify that these 
estimates, reported in Table 1, are consistent with a determination of this 
contribution obtained comparing the lifetimes of decays with the D produced 
in the upstream half of each target with those produced in the downstream 
half of the same target (see Ref. [11] for more details on the target setup). 
Each partition represents a different mixture of hadronic absorptions of decay 
daughters and D mesons. 

The acceptance could be another source of uncertainty. We analyzed this effect 
determining the D° lifetime without the correction function f(t') and removing 
the fiducial volume cut from the set of analysis cuts (this makes the correction 
function almost flat). We obtained results in good agreement with the reported 
values. This check is not possible for the D + because of the longer lifetime, a 
geometric acceptance correction is always needed. A study was performed in 
FOCUS (see Ref. [10] for more details) comparing the acceptance part of the 
Monte Carlo correction with the high statistics Kg — > 7r + 7r~ decays. The result 
of this study showed an excellent agreement between the acceptance observed 
in the data and the acceptance simulated by the Monte Carlo; however we 
assess a 2% uncertainty due to the finite statistics. This 2% uncertainty in the 
f(t') correction function gives a 0.21% and 0.52% uncertainty in the lifetime 
of D° and D + respectively. 

The finite Monte Carlo statistics give a negligible contribution to the system- 
atic uncertainty. 

Table 1 shows the contributions of each of these sources to the total systematic 
uncertainty. For the combined D° lifetime the systematic error (also shown in 
Table 1) is obtained combining the individual sources of systematic uncertainty 
from D° — > K~ir + and D° — > K~it + it~it + . We assume the absolute time 
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scale and the absorption correlated. To obtain the final systematic error the 
uncertainties from the different sources are then added in quadrature. 

The final lifetime values are 409.62±1.15 (statistical) ±1.55 (systematic) fs for 
D° (weighted average) and 1039.42±4.28 (statistical) ±6.97 (systematic) fs for 
D + . A final check was performed for the D° lifetime. We compute the lifetime 
using a combined likelihood, that is forming a global likelihood for K~tt + and 
K~ir + ir~ir + . The fit parameters are the two distinct backgrounds and one 
lifetime. The lifetime from the combined likelihood, 409.62 ±1.15 (statistical), 
is identical to the reported value. 

This measurement of the D° lifetime value is in very good agreement with the 
result we obtained in our lifetime difference paper [4]. 

In conclusion we have measured the lifetimes of the D° and D + mesons. Our 
results are reported in Table 2 along with a comparison with the most recent 
published measurements. 



Table 2 

Measured lifetimes (xlO -12 s). 



Experiment 


D° 


D+ 


E687 [8] 
CLEO II [12] 
E791 [13] 

This measurement 


0.413 ±0.004 ±0.003 
0.4085 ± 0.004llo;oo34 
0.413 ± 0.003 ± 0.004 
0.4096 ± 0.0011 ± 0.0015 


1.048 ±0.015 ±0.011 
1.0336 ± 0.022llo;oi2? 

1.0394 ± 0.0043 ± 0.0070 



Our results will significantly decrease the errors on the current world average 
values for the D° and D + lifetimes. 



From our measurements of the D° and D + lifetimes we can update the de- 
termination of the ratio t(D + )/t(D°): 2.538 ± 0.023. This result and the 
inclusive semileptonic branching ratios [3], D + — > eX = (17.2 ± 1.9)% and 
D° — > eX = (6.75 ± 0.29)%, show that the D° and D + semileptonic decay 
widths are nearly equal: 

T{D+ -> eX) B{D+ -»• eX) r(D°) V ; 



This implies that differences in the total decay widths between D° and D + 
must be due to differences in the hadronic decay sector. 
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Fig. 1. (a) K~tt + invariant mass distribution, (b) K^tt + tt + tt^ invariant mass dis- 
tribution, (c) K~ir + ir + invariant mass distribution. The fits (solid curves) are de- 
scribed in the text and the numbers quoted are the yields. The vertical dashed lines 
indicate the signal region, the vertical dotted lines the sideband. 
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Fig. 2. Decay mode D° — > K~tt + : (a) the correction function /(i'), deviation from a 
flat distribution represents the correction from a pure exponential function; (b) the 
predicted events (histogram) are superimposed to the observed events (points), the 
shaded distribution shows the t' distribution of the background; (c) the background 
subtracted and f(t') corrected t' distribution, the superimposed straight line is an 
exponential with the fitted lifetime. 
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Fig. 3. Decay mode D° — > K~7r + Tr~-K + : (a) the correction function /(<'), deviation 
from a flat distribution represents the correction from a pure exponential func- 
tion; (b) the predicted events (histogram) are superimposed to the observed events 
(points), the shaded distribution shows the t' distribution of the background; (c) 
the background subtracted and f(t') corrected t' distribution, the superimposed 
straight line is an exponential with the fitted lifetime. 
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Fig. 4. Decay mode D + — > K _ 7t + 7t + : (a) the correction function /(£')> deviation 
from a flat distribution represents the correction from a pure exponential func- 
tion; (b) the predicted events (histogram) are superimposed to the observed events 
(points), the shaded distribution shows the t' distribution of the background; (c) 
the background subtracted and f(t') corrected t' distribution, the superimposed 
straight line is an exponential with the fitted lifetime. 
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Fig. 5. The fitted lifetime versus the tjoi detachment cut for: (a) D° — > K~ir + , 
(b) _D° — > K^7r + 7r + 7r~ and (c) D + — > i"C~7r + 7r + . The horizontal dashed lines show 
the interval corresponding to the chosen lifetime ±la. 
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Fig. 6. Invariant mass distributions obtained using tight cuts for: (a) K~ir + 
K~it + it + it~ and (c) K~it + it + . The functions used to fit the data (solid curves 
similar to those of Figure 1 and the numbers quoted are the yields. 
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